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Abstract 

The B — > K*£~£ + (I = fj,, r) is analyzed in a minimally extended Stan- 
dard Model in which the Wilson coefficients have new CP-odd phases. The 
sensitivity of the CP asymmetry and lepton polarization asymmetries to 
the new phases is discussed. It is found that the CP asymmetry is sensitive 
to the new phase in the Wilson coefficient Cj whereas the normal lepton 
polarization asymmetry is sensitive to the phase in the Wilson coefficient 
Cio- Additionally, the correlation between the CP and normal lepton po- 
larization asymmetries is studied. A simultaneous measurement of these 
two asymmetries can be useful in search for the existence of the new sources 
of CP violation beyond the Standard Model. 
PACS: 13.20.He, 11.30.Cp, 12.60.-i 
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I. INTRODUCTION 



Violation of the CP symmetry has now become a well-established fact in Kaon system 
IJ. Once the era of 5-factories start with the operation of KEK-B, B-TeV, LHC-B and 
SLAC's Asymmetric B factory, it will be possible to test the standard model (SM) at one- 
loop accuracy. In general, possible incompatibility of the experimental data with the SM 
predictions will mark the existence of 'new physics' contributions. Among all, an exper- 
imental determination of the CP-violating quantities and their comparison with the SM 
predictions will be particularly useful in search for the new physics effects. 

From the experimental perspective the exclusive decay modes (such as B — > K*^ 0, 
B — > K*£ + £~ , B — > K£ + £~) are easy to measure. From the theoretical view point, however, 
the corresponding inclusive modes (6 — > S7 and b — *> s£ + £~) can be cleanly estimated as the 
only machinery needed are the Wilson coefficients describing the short-distance physics. A 
proper description of the exclusive decay modes, on the other hand, depends on both Wilson 
coefficients (short-distance physics) and the hadronic form factors (long-distance physics). 
This causes a relative increase of the uncertainties due to hadronization effects. 

For the purpose of studying the sources of CP violation, it is convenient to concentrate on 
those observables which are sensitive to the possible CP phases. Among these, for example, 
CP asymmetries and lepton polarization asymmetries are such ones ||- |J. Recently, a 
detailed study of the lepton polarization asymmetries in B — > X s £ + £r decay has been 
performed in a rather general model by including nine additional Wilson coefficients not 
found in the SM ||. Keeping this kind of short-distance generality it is convenient to 
discuss the exclusive decay modes such as B — > K*£ + £~ J7|||. Such an analysis will be 
useful for a first-hand comparison with the experiment as the inclusive modes are generally 
hard to measure. 

In Sec. II we start with a general non-standard description of the short-distance physics 
as in |J. Then we parametrize the long-distance quantities by appropriate form factors and 
obtain the hadronic transition amplitude. In Sec. Ill we derive general analytic expressions 
for asymmetries will be give. In doing this all sources of CP violation will be ascribed 
to short-distance physics. In Sec. IV the asymmetries and their relation to the Wilson 
coefficients will be analyzed numerically. In Sec. V results are discussed and the conclusion 
is stated. 
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II. THE DECAY AMPLITUDE 



M(6 -> s£ + £~) = -KzV th Vl{ C SL sia^(m s L)b£^£ + C BR sia ^{m h R)b £7 J 



The exclusive B decays B — > K*£ + £~ are conveniently described by the partonic decay 
b — > s£ + £~ at distances (^(M^ 1 ). The effective Hamiltonian describing this rare b decay at 
the scale \i ~ My/ should, however, be evolved down to mesonic mass scale fi ~ mj using 
the QCD evolution equations. Then the decay amplitude describing b — > s£ + £~ takes the 
form H[7] 

-T^^tfeKi Csl sia^ — (m s L)b £7/ + C BR sia^ — 
V27r ( <r <r 

+ Crr sja^R £r^£r + Sifoij £ L £ R + CijLLi? s R b L £ L £ R 

+ CWl s L 6 K I R £ L + C RLRL s R b L £ R £ L + C T sa^bla^t 

+ iC TE e^sa liV b£a aP £^ , (1) 

where each of the Wilson coefficients C$l, • • ■ , Cte is evaluated at the 5-meson mass scale, 
ji ~ m b . In this expression, L(R) = (1 — (+)7s)/2 are the left (right) projection operators, 
are the elements of the CKM matrix, and q = p B —px* — P+ +P- is the momentum transfer 
to the dilepton channel. This decay amplitude has a rather general form as it includes 
nine additional operators not found in the minimal standard model. The only simplifying 
assumption about this decay amplitude will be twofold: (1) Neglect of the strange quark 
mass everywhere in the analysis, (2) Neglect of the tensor operators having the coefficients 
Ct and Cte- The former is justified by the smallness of the ratio m s /mb and the latter 
is justified by the previous analyzes which show that their contributions are much smaller 
than other operators (for details, see [[|). 

The quark level decay amplitude ([I]) controls the semileptonic decays B — > (K, K*) £ + £~ . 
The amplitudes for these exclusive decays can be obtained after evaluating the ma- 
trix elements of the quark operators in ([!]) between the \B{p B )) and (K*(pK*)\ states. 
In particular, explicit expressions for (K* |s7^(l ± 75)6] B), {K*\sia llv q v {l+^)b\B) and 
(K* \s(l ± 75)6! B) are needed. Computation of such hadronic matrix elements is bound to 
parametrizations of the form factors depending only on the momentum transfer square, or 
equivalently, the dilepton invariant mass m 2 u = (p B — Pk*) 2 = (p+ + P-) 2 = q 2 ■ Introducing 
appropriate from factors one obtains 



(K*(p K *,e)\s llM (l± l5 )b\B(p B )) = 

p p K «q a ± Kims + m K ,)A 1 (q 2 ) T ^{PB + p K *)»(e*q) ^ 
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where the explicit expressions for V(q 2 ), A ^ 2t3 (q 2 ) and Tx^{q 2 ) will be given below. 

To ensure the finiteness of (0) as q 2 — ► 0, it is usually assumed that A 3 (q 2 = 0) = A (q 2 = 
0). Besides, to calculate the matrix elements of the scalar operators, (K* \s(l ± 75)6] -B), it 
is necessary to contract (|2|) with q^ and use the equation of motion, giving 



(K*(p K .,e) \s{l ± 75 )6| S(p B )) = — { T i(s*q)(m B + m K ^A 1 {q' 2 ) 

± z(m B - m x *)(£*g)A 2 (? 2 ) ± 2im K *(s*q) U 3 (g 2 ) - A (q 2 ] 



}■ 



(4) 



Additionally, again using the equation of motion, the form factor A 3 can be expressed as a 
linear combination of the form factors A\ and A 2 (see [H]]) 
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7ns - thk* 



2rriK* 2rriK* 
Having this relation at hand, one finally obtains 



{K*(p K .,e) \s(l±y s )b\B(p B )) 



— ( =1= 2im 



*q)A (q 2 )} • 



(5) 



(6) 



This completes the evaluation of the necessary transition matrix elements. 

As mentioned before the form factors entering (@)-(f]) represent the hadronization process 
which lacks a Lagrangian description. They are thus generally computed in framework of 
certain nonperturbative approaches such as chiral theory f[T| |, three point QCD sum rules 



method ||10|| , relativistic quark model by the light-front formalism |T!|, effective heavy quark 
theory Jl3j and light cone QCD sum rules fl4Hl6|. In what follows we will use the results 
of the work in which the form factors are described by a three-parameter fit where the 
radiative corrections up to leading twist contribution and SU(3)-breaking effects are taken 
into account. Letting F(q 2 ) e {V (q 2 ) , A (q 2 ) , A^q 2 ) , A 2 (q 2 ) , A 3 (q 2 ) , T^q 2 ) , T 2 (q 2 ) , T 3 (g 2 )} , 
the g 2 -dependence of any of these form factors could be parametrized as 
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At 


-+K* 


0.34 ±0.05 
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-0.023 


Af 


->K* 


0.28 ±0.04 


1.18 


0.281 


yB- 


■+K* 


0.46 ±0.07 
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0.575 


T B- 


+K* 


0.19 ±0.03 


1.59 


0.615 


rpB- 

1 2 


■>K* 


0.19 ±0.03 


0.49 


-0.241 


rpB- 

J 3 


■>K* 


0.13 ±0.02 


1.20 


0.098 



TABLE I. The form factors for B — ► K*£ + £ in a three-parameter fit. 

where the parameters F(0), ap and bp are listed in Table 1 for each form factor. Here 
s = q 2 /m 2 B is the dilepton invariant mass in units of S-meson mass (See, |15|,16]|). 

Making use of the hadronic matrix elements of the basic quark current structures 

in it is straightforward to determine the decay amplitude for B — > K*£ + £~ decay: 

Ga 



M(B -> K*l 



4^'* 

x - 7sK [ - 2V Ll e flupa e* u p p K ,q CT - %V h2 e^ + iV L3 (£*q)(p B +Pk*)» + iV LA (e*q) qil 

+l lfl (l + 7sK [ - 2V Rl e^ pa e* u p p K «q a - iV^e* + iV R3 (s*q)(p B + p K *)^ + iV Ri (e*q)q p 
+1(1 - l5 )£\iS L (e*q)] + 1(1 + l5 )£\iS R (e*q) 



(7) 



where Vi^ and V Ri are the coefficients of left- and right-handed leptonic currents with 
vector structure, respectively. Clearly, Sl,r are the weights of scalar leptonic currents with 
respective chirality. These new coefficients are functions of the Wilson coefficients in the 
partonic decay amplitude (|I|) and the form factors introduced in defining the hadronic matrix 
elements above. Their explicit expressions are given by 
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where q 2 dependencies are implied. It is clear that all these effective form factors are 
functions of the specific form factors and the Wilson coefficients in ([!]). Therefore, 

they carry information about both short- and long-distance physics. 

The hadronic matrix element M(B — ► K*£ + £~) is the basic machinery for the computa- 
tion of all physical quantities pertaining to this decay. In particular, the computation of the 
energetic distributions, total rate and various asymmetries follow from (B — > K*£ + £~) using 
the usual methods. In the next section, necessary asymmetries and other relevant quantities 
will be computed. 



III. ASYMMETRIES 



For an analysis of the asymmetries it is necessary to compute the differential decay rate 
for B — > K*£ + £~ decay. For unpolarized leptons at the final state, using the decay amplitude 
in (|7|), the differential decay rate is found to be 

/ dT \ G 2 a 2 



V tb V*\ 2 X 1 ^ 



dq 2 J Q 2 14 vr 5 m B 

-h m %s - m 2 )(\V Ll \ 2 + \V Rl \ 2 ) + 2m, 2 Re(V Ll V^; 



x I 32Am 4 
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+ 96m, 2 Re(V L2 V^ 2 ) - -m 2 B m e X Re[(V La - Vr 2 )(S* l - S* R )] 

+ ^m|m 2 ARe[V! 2 (V L4 + V R[i - V Ri ) + V* R2 (V Ls - V Li + V Ri ) - (S L S* R ) 

+ -m Bmi {l - r)X {Re[(V L3 - V R ,)(S* L - S R )}} 

+ -m B m\{\ - r)X {Re[-(V L3 - Vr 3 )(V* L4 - V R J]} 

H m B m 2 X(2 + 2r - s) Re(Vz 3 V£ 3 ) - -m B m iS X Re[(V L4 - Vr 4 )(S* l - S* R )\ 

2 



r 

4 4 2 x 

-m R mi,sX 
r 



\VlS + \VrX - 2Re(V L4 V^)j + -m B (m B - 2m\)X [ \S L \ Z + \S R \ 
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3rs 
4 
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rs 



m 2 B \ 



m\{2 — 2r + s) + m 2 B s{l — r — s) 



2m\{\ - 6rs) + m B s(X + 12rs) 



Srs 



B X [m%s\ + m] [2A + 3s(2 + 2r 



_ Re(V Z2 V2 3 )+Re(V ii2 V^) 
|V i2 | 2 +|V % | 

)}}[\Vl 3 \ 2 + \V R3 \ 



(8) 



where the subscript "0" is intended for the unpolarized decay rate. In this expression 
s = q 2 /m 2 B , r = m 2 K ,/m 2 B , v 2 = l — (4mf)/q 2 , and finally A(l, r, s) = l + s 2 + r 2 — 2r — 2s — 2rs 
is the familiar triangle function. 

Our next task is the calculation of the lepton polarization asymmetries. For the computa- 
tion of these asymmetries the unpolarized decay rate (|8]) is not sufficient. The measurement 
of these asymmetries require the specification of the total number of leptons of a given kind 
(for example, negatively-charged) in a given direction. Therefore, it is necessary to take into 
account the polarization of the lepton beam in a given direction. Considering, for example, 
the negatively- charged lepton, one can introduce the following three polarization vectors in 
the rest frame of l~: 
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(9) 



where efj ■ e} 



Si j, P- 



\(Pk* x p- 

ti = 0, i, j = L, T, N. Here, e^, and e/v correspond, respectively, to 
the 'longitudinal', 'transversal' and 'normal' polarization directions of l~ with respect to its 
direction of motion, e L . One notices that and are co-planar, and e/v is perpendicular 
to this plane. In the rest frame of £~ , the temporal components of the corresponding 
four-vectors vanish. However, in the dilepton rest frame (that is, q = 0), the four-vector 
corresponding to e*x is boosted to (pl/m^, {Eijm^ei) leaving e T and unchanged. In the 
following, all results will be conveniently given in the dilepton rest frame. 

The differential decay rate for any spin direction n of the £~, where n is a unit vector in 
the l~ rest frame satisfying n - n — l,n ■ p~ — 0, can be expressed in the following form 

dT{n) 1 ( dV 



dq 2 



dq 2 



1 + (P L e L + P N e N + P T e T 



■ n 



(10) 



where the coefficients of unit vectors, Pl, Pn and Pt, are recognized as the 'longitudinal', 
'normal' and 'transversal' polarization asymmetries. A simple formula for extracting these 
asymmetries from the polarized decay rate follows from ( |10D itself: 
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where % — L,T, N. One notes that the denominator in this expression is identical to the 
unpolarized decay rate (||). On the other hand, the numerator depends on the spin direction 
of the lepton under consideration. In essence what Pi{q 2 ) measures is the difference between 
the rates for a particular direction and its opposite for a given dilepton invariant mass 

^2 



me = VQ 

Using the hadronic decay amplitude (^) in the general polarization asymmetry formulae 
(p!l|), after a lengthy calculation the polarization asymmetries Pl, Pn and Pt are found to 
have the following explicit expressions: 
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(12) 



where A is the expression within the curly parenthesis in the unpolarized differential decay 
rate in These expressions for the polarization asymmetries are quite general except for 
the neglect of strange quark mass and the tensor operators (the last two operators in ([!])) 
as mentioned before. 

Before proceeding, it is convenient to make a few useful observations on the lepton 
asymmetries. Particularly interesting one is the massless (light) lepton limit: m, — > 0. In 
this case, Pl depends only on the bilinears of and V Ri , that is, the effects of the scalar 
operators in (p]) completely decouple. On the other hand, Pt and Pn happen to depend only 
on the interference terms between the coefficients of the vector operators (V^ and VrJ and 
those of the scalar operators (Sl,r)- However, one notices that the scalar operators in (|I|) can 
be induced by an exchange of the scalar particle (such as two Higgs doublet models |I7| ) in 
which case the coefficients Sl,r are necessarily proportional to the lepton mass. Therefore, 
in the limit of massless (light) leptons only the longitudinal asymmetry Pi can remain non- 
vanishing. Conversely, in near future, if experiment yields non-vanishing Pjy and Pt for 
B — > K*e + e~ decay this would imply the generation of scalar operators by mechanisms 
beyond the Higgs model where the fermion scalar coupling is always proportional to the 
fermion mass. 

In general, the lepton polarization asymmetries ( |T2"D are able to probe real as well as 
imaginary parts of the effective form factors V^, Vr^ and Sl.r. As the parametrization 
(f|]) shows the hadronic form factors are inherently real, and thus the imaginary parts of 
Vij> Vr 4 and Sl,r in (|7|) can come only from the Wilson coefficients in ([I]). Below we 
will keep this picture, that is, we will be dealing only with the CP violation effects due to 
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short-distance physics parametrized by the Wilson coefficients. At this point it is useful 
to distinguish between CP properties of the phases in the Wilson coefficients. In principle, 
Cbr, ■ ■ -Cte all can have finite phases; however, these phases can have strong and weak 
subparts. Here by strong and weak we mean even and odd phases under CP conjugation. To 
be able to distinguish such distinct components of the phases it is not sufficient to analyze the 
polarization asymmetries alone. One, in particular, has to consider the CP asymmetry of the 
decay which is inherently sensitive to CP character of the phases of the Wilson coefficients. 
Using the unpolarized decay rate (§), the CP asymmetry for B — > K*£ + £~ decay is defined 
by: 



A CP (q 2 ) 






(13) 



where the processes to which dT/dq 2 refers are explicitly shown in the arguments. Making 
use of the explicit expression for the unpolarized decay rate (|8|) one can determine the 
detailed dependence of Acp(q 2 ) on the model parameters. For this purpose it is useful to 
introduce the following parametrization for the quantities V^, and Sl,r (0): 
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where i — 1, ■ ■ -,4. In this expression subscript "s" ("to") stands for strong (weak) phases 
mentioned above. By definition, V's and 5's are combinations of hadronic form factors 
and Wilson coefficients so that the phases <p W)S defined by (|14]) are explicit functions of the 
dilepton invariant mass. With this definition of the from factors it is possible to find a 
suggestive form for the CP asymmetry: 
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where A0"' b = 0™ — 0^. The quantity £ in the denominator is even under all these phases, 
and has the expression 
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Until now the decay rate (Bfi, the lepton polarization asymmetries (O) and CP asym- 
metry ([13j) have been computed by adopting a rather general quark level amplitude ([!]) for 
I? — ► K*£ + £~ decay. Presently this exclusive decay has a direct bound coming from recent 
CDF measurement |TB[: BR(B — > K*fi + fi~) < 4.0 x 10~ 6 . In addition to this direct bound, 
existing CLEO result for BR(5 — > X*7) imposes another important, albeit partial, con- 
straint on the parameter space. Indeed, using the notation of ([J) and appropriate form 
factors derived in (§)-© the total decay rate for B —>■ K*^ can be written as 



Y{B -> K*i) 



G 2 am\ml 



1287T 



""""vu^fi-^lVl^ 



(17) 
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which constrains directly \Cbr Ti(0)|. Therefore, the norm of the Wilson coefficient for the 
dipole operator is fixed by BR(i? — > K*^). This constraint, however, does not say anything 
about the CP violation potential of Cbr- 

IV. LARGE CP PHASES WITHIN AN SM LIKE MODEL 

As as been emphasized before, the underlying model for the partonic decay amplitude 
(|I|) is rather general one. In principle, all the Wilson coefficients can be non-zero and may 
have arbitrary phases. However, a realistic model should meet the requirements of the SM 
to leading order, and especially, should not cause unacceptable deviations from the existing 
experimental data confirmed already by the SM. For this reason it is convenient to establish 
the discussion on a general model with close reference to the SM predictions. 

Therefore, here we follow a minimal prescription such that the general Wilson coefficients 
in ([!]) are 

• endowed with new phases beyond the SM, 

• identical to SM ones when these phases vanish. 

Such an approach obviously neglects the new physics contributions to the norms of the 
Wilson coefficients; however, at the aim of determining the information content of the po- 
larization and CP asymmetries on the sources of CP violation, it suffices. This is kind of a 
minimal approach for parametrizing the new physics CP violation for asymmetry measure- 
ments in B — > K*£ + £~ decay. 

Adopting this approach, one can make the following assignments for the general Wilson 
coefficients in (JiJ). First, the coefficients describing the scalar-scalar type interactions vanish 
identically 

Clrrl = Crllr = Clrlr = Crlrl = 0. (18) 

It is known that such coefficients exist, for example, in the two Higgs doublet models (2HDM) 
which have an extended Higgs sector compared to the SM. In such models these scalar-scalar 
interactions are induced by the Higgs exchange, and the resulting Wilson coefficients are 
proportional to m&m^/m^ which is maximal for £ — r. However, to the extent one neglects 
m s /m&, mbme/rri^ is, too, negligible in the light of recent LEP limit on the Higgs boson 
mass rrih |BJ is taken into account. 

The Wilson coefficient for the dipole operator Cbr obeys 
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C BR = -2C e 7 ff (m b ) = -2 (c 7 (m b ) - ^C 5 (m b ) - C 6 (m b )^j (19) 

where 07 is an arbitrary phase, and it is not constrained by BR(5 — > K*^f) at all. 
Finally the coefficients of the vector-vector interactions in ([!]) are given by 

C LL(LR) = C e 9 ff (m b ) e 1 ^ - (+)C 10 (m b ) , C RL = C RR = (20) 

where the coefficient Cio is known to be scale independent: Cio(m b ) = Cio(Mjy). 

In the SM the Wilson coefficients (m b ) and Cio(m b ) are strictly real as can be read 
off from Table II. Moreover, the SM prediction for Cj (m b ) is already consistent with the 
CLEO determination of BR(5 — > X*7) in flT7|). Therefore, through the choice of C R r 
in (|19D the experimental constraint is already taken into account. Although individual 
Wilson coefficients at \i ~ m b level are all real (see Table II) the effective Wilson coefficient 
C 9 \m b ,q 2 ) has a finite phase, and is an explicit function of the dilepton invariant mass, 
q 2 . To see its phase content it is useful to reproduce its explicit expression here: 

C e 9 ff (m b ) = C 9 (m b ) jl + (s) j + Y SD (m b , s) + Y LD (m b , s) (21) 

where C 9 (m b ) is read off from Table II. Here u (s) represents the 0(a s ) corrections coming 



from one-gluon exchange in the matrix element of the corresponding operator 9 |20 



oo (s) = --TT 2 - -Li 2 (s) - - In (s) In (1 - s) 5 + 4 \ In (1 - s) 

K 1 9 3 V ; 3 V ; V ; 3(1 + 2s) V ; 

_ 23(1 + 3) (1-23) S + 93-63 2 

3(l-s) 2 (l + 2s) 1 ; 3(l-s)(l + 2s) ' 1 ; 

In (^1|) Ygx) and Yld represent, respectively, the short- and long-distance contributions 



of the four-quark operators Oi=i ) ...,6 p0| , pi| . Here Ysd can be obtained by a perturbative 
calculation 

Isd K, s) = 5 (™o s) [3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ] 
-^(1,S)[4C 3 + 4C4 + 3C 5 + C 6 ] 

- l -g (0, 5) [C 3 + 3C 4 ] + [3C 3 + + 3C 5 + C 6 ] 

- [ 3Cl + ^ b (°» 5 ) - 5 (™c, *)] > ( 23 ) 

where the loop function g (m q , s) stands for the loops of quarks with mass m q at the dilepton 
invariant mass s. This function develops absorbtive parts for dilepton energies s = 4m 2 : 
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^ 8 > * 8 4 2 / \ a r 

g (m q , s) = --\nm q + — + -y q - - (2 + y q ) - y q \ 

x{9(l-y fl ) (in - + y/ IHi - m) + 6(y g - 1)2 arctan— L= }, (24) 

V 1 - V 1 - Vq J VVq ~ 1 J 

where ?n 9 = m q /m b and y g = 4rhg/s. Hence, due to the absorbtive parts of g (m q , s), there 
are strong phases coming from Y$d- One, in particular, notices the terms proportional to 
g (0, s) which have a non-vanishing imaginary parts independent of the dilepton invariant 
mass. 

In addition to these perturbative contributions, the cc loops can excite low-lying char- 
monium states ip(ls), ■ • • , ip(6s) whose contributions are represented by Y LD |22]| : 

x •£ ™*r(v^ e n-)M v (25) 

v^{is),-,^ s ) (M$. - sm 2 b - iM v V Vi ) 

where the Fudge factors for B -»• K*Vi -> transition, and C(°> = 3Ci + C 2 + 

3C3 + C4 + 3C5 + Cq. Here the sum runs over all all charmonium resonances with mass my i 
and total decay rate T^. Contrary to Y$d, the long-distance contribution Yld has both 
weak and strong phases. The weak phases follow from the CKM elements whereas the strong 
phases come from the s values for which i-th charmonium state is on shell. Therefore, the 
Wilson coefficient Cg^ (m&) has both weak and strong phases already in the SM. 



Ci(m&) 


C 2 {m b ) 


C 3 (m b ) 


C 4 (m b ) 


C b (m b ) 


Ce(m b ) 


C e 7 ff (m b ) 


C 9 (m b ) 


C w (m b ) 


-0.248 


1.107 


0.011 


-0.026 


0.007 


-0.031 


-0.313 


4.344 


-4.669 



TABLE II. The numerical values of the Wilson coefficients at /i ~ m b scale within the SM. The 
corresponding numerical value of C° is 0.362. 



In this sense, the Wilson coefficients C% (m b ) and Cio(m b ) cannot develop any strong 
phase, and thus, 07 and 0io should necessarily originate from physics beyond the SM. A 
few observations on the asymmetries help much in simplifying the analysis: (i) Due to the 
dependencies of the asymmetries on the Wilson coefficients, it is clear that one can re-phase 
one of the Wilson coefficients. For instance, one can choose 0g = leaving Cg^(m b ) with 
its SM phases only, (ii) As mentioned above, the Wilson coefficients (m b ) and Cio(m b ) 
cannot develop strong phases from light quark loops so that 07 and 0io can be chosen to 
have purely weak character. 
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In the light of analytic derivations as well as particular observations mentioned above, 
one can investigate the dependence of the asymmetries on these new weak phases 07 and 

010 to have an estimate of their information content. It is conceivable that such an analysis 
will provide a tool to mark possible sources of CP violation beyond the CKM matrix. 

V. NUMERICAL ESTIMATES 

In this section we present our numerical estimates for the asymmetries A C p, Pl, Pt and 
P N for B — > K*fi + fi~ and B — > K*t + t~ decays separately. We take hadronic form factors 
from Table I and the Wilson coefficients from Table II. For the remaining parameters we 
take m b = 4.8 GeV, m c = 1.35 GeV, m B = 5.28 GeV, m K « = 0.892 GeV. 

The dilepton invariant mass has the kinematical interval Am] < q 2 < {m B — m K *) 2 in 
which the charmonium resonances can be excited. The dominant contribution comes from 
the three low-lying resonances J '/ip,ip' \ip" in the interval 8 GeV 2 < q 2 < 14.5 GeV 2 . In 
order to minimize the hadronic uncertainties we will discard this subinterval in the analysis 
below by dividing the q 2 region to low and high dilepton mass intervals 

Region I : Am] < q 2 < 8 GeV 2 , 

Region II : 14.5 GeV 2 < q 2 < (m B - m K ,f, (26) 

where the contribution of the higher resonances do still exists in the second region. 

Due to 1/q 2 factor in front of C BR , in Region I the contribution of the dipole type 
operator dominates. Therefore, asymmetries which involve the differences of the decay rates 
are suppressed in Region I compared to ones in Region II. This property will be illustrated 
in Figs. 1-2 and the remaining analysis for the asymmetries will be performed only for 
Region II where the asymmetries are larger. 

As mentioned previously, in the model under concern, there are two weak phases: 07 
and 0io- However, a close inspection of the CP asymmetry shows that, it is independent of 
0io- This follows from the fact that CP asymmetry can exist only when interference terms 
involve strong and weak phases. In this model, similar to SM, there is no interference terms 
involving Cg^ and C w . For this reason A C p is independent of 0i O . 

First, we illustrate A C p in fa-q 2 plane for B — > K*fi + fi~ decay for Region I and Region 

11 in Fig. 1 and Fig. 2, respectively. In both figures we take 0i O = 0, and as noted above 
Acp is already independent of this phase. In Region I the CP asymmetry is practically 
independent of q 2 , and becomes maximal for marginal CP violation, 7 = tt/2. In Region II, 
however, the q 2 dependence is comparatively enhanced as the dominance of dipole coefficient 
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is now reduced. Besides, as figures suggest the CP asymmetry in Region II is one order of 
magnitude larger than in Region I, and this confirms our expectation above. 

One notes that the average asymmetries could be measured more easily in experiments. 
Therefore, from now on we will discuss averaged CP and lepton polarization asymmetries 
in Region II. The averaging procedure is defined by 

(Q) = Jl4 ; 5GeV2 ^ dq2 q (27) 

J 14.5 GeV 2 dq 2Ui l 

where Q = P L , P N , P T or A C p. 

Depicted in Fig. 3 (Fig. 4) is the 07 dependence of the averaged asymmetries (Pl), (Pr), 
(Pn) and (Aqp) at 1O = ( 0io — tt/2) for B — > K*fi + fx~ decay. Similarly, Figs. 5 and 6 
show the 07 dependence of the same quantities for B — ► K*t + t~ . As noted before, the CP 
asymmetry depends only on 07; however, as these figures show clearly among all asymmetries 
Pjv is very sensitive to 4>iq'- For 0i O = (0io — tt/2) Pn is purely positive (negative). In 
addition to this, Pn at 0io = tt/2 is one order of magnitude larger than that at 0io = 0. 
This property is valid for both and t + t~ final states. Besides, since Pn is proportional 

to the lepton mass, the B — > K*t + t~ decay is much more relevant for its measurement. 
This sensitivity of Pn on 0i O can be explained as follows: Pn depends on the imaginary 
part of the bilinear combinations of the Wilson coefficients, such as Im 
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When 

0io = tt/2 ( 0io = 0) Cio is pure imaginary (real) and therefore Im[Cg /J r C* ] = |Ci |Re[Cg^] 
{\m[Cl ff C* w } = |Cio|Im[C 9 e// ]). Since |Re[C 9 e// ]| » |Im[C 9 e// ]|, P N at 1O = tt/2 is roughly 
one order of magnitude larger than its value at 0io = 0. Remaining two asymmetries, Pl 
and Pt, are less sensitive to 0io- 

In Fig. 7 and 8 we present the correlation between {Aqp) and (Pn) for B — > K*t + t~ 
decay by varying 7 from to 2tt at 0i O = and 0io = tt/2, respectively. For channel 
Pn is much smaller so we do not analyze this case. The SM predictions are given by the 
intersections of (A C p) = line and the curves themselves. Due to the sign ambiguity of C 7 
there are two solutions. All other points on the curves are generated by the new physics 
phases. If a simultaneous measurement of (Acp) and (Pn) gives a point on the curve and 
if this point is distinct from the SM prediction then this will be an indication of the new 
physics contribution. Moreover, such a simultaneous measurement enables us to determine 
the sign of the new phases unambiguously. 
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VI. CONCLUSION 



In this work we have adopted a model-independent approach in studying the sensitivity 
of the CP and lepton polarization asymmetries to new CP phases. In parcticular, we have 
taken the Wilson coefficients being identical to the SM ones except for their phases. The 
main result of the present study is that the CP asymmetry and normal lepton polarization 
asymmetry are the most sensitive quantities to new sources of weak phases beyond the SM. 
While A C p is sensitive to 07 only, P N is more sensitive to Therefore, measurement 
of these two asymmetries can establish the existence or absence of the new sources of CP 
violation beyond the SM. Moreover, a simultaneous measurement of the averaged CP and 
normal polarization asymmetries will unambiguously determine the sign of the new phases. 

In a specific model such as 2HDM or supersymmetry the Wilson coefficients possess 
new CP phases not found in the SM. The question of how infromative the asymmetries in 
B — > K*l + £~ decay about new sources of CP violation in 2HDM model or supersymmetry 
will be discussed elsewhere. 
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FIGURE CAPTIONS 

Fig. 1 The dependence of the CP asymmetry A C p for B — > K*fi~fi + on q 2 and 7 at 
0io = for Region I. 

Fig. 2 The same as in Fig. 1 but for Region II. 

Fig. 3 The dependence of (A CP ), (P L ), (P T ) and (P N ) for B -> on 7 at <f> 10 = 

for Region II. 

Fig. 4 The same as in Fig. 3 but for 1O = tc/2. 

Fig. 5 The same as in Fig. 3 but for B — > K*t~t + decay. 

Fig. 6 The same as in Fig. 5 but for 1O = ir/2. 

Fig. 7 The correlation between the averaged CP and normal lepton polarization asym- 
metry at 0io = for B — > K*t~t + decay. 

Fig. 8 The same as Fig. 7 but for 1O = ir/2. 
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FIG. 1. 
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